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ABSTRACT 

The effect of galaxy interactions on star formation has been investigated using Data Re- 
lease 1 of the Sloan Digital Sky Survey (SDSS). Both the imaging and spectroscopy data 
products have been used to construct a catalogue of nearest companions to a volume lim- 
ited (0.03 < z < 0.1) sample of galaxies drawn from the main galaxy sample of SDSS. Of 
the 13973 galaxies in the volume limited sample, we have identified 12492 systems with 
companions at projected separations less than 300 kpc. Star-formation rates for the volume- 
limited sample have been calculated from extinction and aperture corrected Ha luminosities 
and, where available, IRAS data. Specific star formation rates were calculated by estimat- 
ing galaxy masses from z-band luminosities, and r-band concentration indices were used 
as an indicator of morphological class. The mean specific star-formation rate is significantly 
enhanced for projected separations less than 30 kpc. For late-type galaxies the correlation ex- 
tends out to projected separations of 300 kpc and is most pronounced in actively star-forming 
systems. The specific star-formation rate is observed to decrease with increasing recessional 
velocity difference, but the magnitude of this effect is small compared to that associated with 
the projected separation. We also observe a tight relationship between the concentration in- 
dex and pair separation; the mean concentration index is largest for pairs with separations 
of approximately 75 kpc and declines rapidly for separations smaller than this. This is in- 
terpreted as being due to the presence of tidally-triggered nuclear starbursts in close pairs. 
Further, we find no dependence of star formation enhancement on the morphological type or 
mass of the companion galaxy. 
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1 INTRODUCTION 

It is well established, both observationally and theoretically, that 
galaxies do not evolve in isolation from one another but that 
their present day structure is the result of sequential merg- 
ers and encounters of ga laxies and proto-galaxies. For example, 
iToomre & Toomre 1 1972) used numerical simulations to show that 
many of the features of peculiar galaxies can be explained by re- 
cent galaxy-galaxy interactions and mergers. Indeed, in the cur- 
rently widely favoured hierarchical models of galaxy formation, 
all of the galaxies at the present e poch are the results of mergers 
and accretion (e.g.. Cole et al. 2000). 

The problem of predicting the outcomes of galaxy encounters 
is difficult. Much progress has been m ade using numerical simu- 
lations (e.g., Mihos & Hemauist '199^, following the pioneering 
work of Toomre & Toomre ( 1972), but several outstanding prob- 
lems remain. The fate of the inter-stellar medium (ISM) and even- 
tual conversion to new stars is particularly difficult to follow due 
to the finite resolution of simulations and the complicated physics 
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of the ISM. In this paper we investigate one aspect of this problem, 
the effect of tidal interaction on the star formation rate of galaxies, 
using data from the Sloan Digital Sky Survey (SDSS: lYork et all 
2000). 

The connection between star formation and galaxy interac- 
tions has been suspected ever since the first large surveys, biased 
towards star forming systems were conducted: first in the ultravi- 
olet continuum, then in emission lines and the far infrared. The 
far infrared survey by the Infrared Astronomical Satellite (IRAS), 
in particular, found a population of star forming systems in the 
nearby universe with LpiR > 1 x lO^'L©, corresponding to an in- 
ferred star formation rate greater than 22M0year~'. These star- 
forming galaxies exhibit, on average, more peculiar morpholo- 
gies and are more likely to be mergers. Furthermore, these ob- 
servational indications of interaction/merger become more com- 
mon and pronounced in systems with higher star-formation rates. 
The most luminous of the IRAS sources, i.e., those with LpjR ^ 
3 X lO'^Lp), being universally classified as strong interactions or 
mergers (.Sanders & Mirabel 1996). 

A number of existing studies have used samples of inter- 
acting galaxies to examine statistically the causal relationship 
between galaxy interactions and star formation. For example. 
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iLarson & TinslevI il978h found a much higher dispersion in the 
U —B versus B ~V colour-colour plot for galaxies selected from 
Arp's atlas (1966) as compared to a control sample, indicating 
an enhanced star form ation rate within the last ~ 10^ years. 
iKennicutt et alj il987h used Ha line and far infrared observa- 
tions to examine the influence of interactions on the global star 
formation of a complete sample of close pairs, and a sample of 
Arp systems. They found a higher than average star formation 
rate in interacting systems, but also found that a large fraction of 
galaxi es in close pair s exhibit normal star formation rates. Sim- 
ilarly, [Sushous^^] 0,988) examined the far infrared emission 
in a sample of strongly interacting galaxy pairs, finding enhanced 
emission in some, but not all systems. Bergvall et al. (2003) have 
looked at star formation in a sample of 59 interacting and merging 
systems, and 38 isolated galaxies, using spectroscopic and pho- 
tometric observations in the optical/near-infrared. In contrast to 
other results, they find that the global UBV colours do not sup- 
port significant enhancement of the star-formation rate in interact- 
ing/merging galaxies. 

Defining samples of interacting or merging galaxies to inves- 
tigate the effects of interaction on star formati on is not straight 
forwar d. One technique, based on the work of iToomre & Toomrd 
is to select galaxies with peculiar morphologies, espe- 
cially those showing tidal tails or galactic bridges. Such s ystem s 
are often drawn from Arp's Atlas of Peculiar Galaxies (19661. 
This technique has been used by a nu mber of author s, including, 
Larson & TinslevI jl97d) and Kennicu tteTai] jl987h . The draw- 
back, as illustrated bv iToomre & Toomr3ti972h themselves, is 
that the degree of morphological disturbance is sensitive to the or- 
bital parameters of the interaction, and so can not be relied upon to 
trace all interactions. Additionally, the features characteristic of in- 
teraction are often of low surface brightness, requiring deep imag- 
ing for detection, and must be identified by visual inspection. 

When the interacting systems are separated by distances com- 
parable to, or larger than, their optical extents, it is relatively easy 
to resolve the galaxies in imaging data and, based on their pro- 
jected separation, identify them as a physical pair. If spectral data 
are available, the velocity separation can also be used. Using this 
technique, survey data can by analysed for close pairs yielding 
large samples of interacting galaxies, albeit usually missing the 
galaxies in the later stages of interaction and merger. 

Large area redshift surveys have ena bled studies of large sam- 
ples of interacting galaxies. For example. lBarton et alJi200ol) have 
analysed optical spectra from a sample of 502 galaxies in close 
pairs and N-tuples from the Cf A2 redshift survey finding the equiv- 
alent widths of Ha anti-correlate strongly with pair spatial and ve- 
locity separation. Similarly, Lambas et al. (2003) examined 1853 
pairs in the 100k public release of the 2dF galaxy survey and find 
star formation in galaxy pairs to be significantly enhanced over 
that of isolated galaxies for separations less than 36 kpc and veloc- 
ity differences less than 100km s^'. 

In this paper we follow a similar approach and use the SDSS 
to define a large sample of 12861 galaxies with identified com- 
panions. Our aim is to establish the relative importance of galaxy 
interactions in determining star formation. 

A number of studies have established a link between den- 
sity of environment and star-formation rate, for example using 
the SDSS (Gomez et al. 2003), 2dF (Lewis et al. 2002), and Las 
Campanas Redshift Survey (Hashimot o et alj[l998 ). In all three 
studies a correlation is found between density of environment and 
star-formation rate in the sense that o ngoing star formatio n is 
suppressed in the denser environments. iGomez et alJ i2003tl and 
iHashimoto et aL (.1998.) have attempted to decouple this effect 



from the densi ty-morpholog y J|Dres sleJll98(l) or radius morphol- 
ogy relations jWhitmore etalll993h finding that the density/star- 
formation rate relationship exists independent of the density- 
morphology relationship. 

The structure of this paper is as follows: In Section 2 we dis- 
cuss the compilation of a volume- and luminosity-limited sample 
drawn from the SDSS, the method used to define galaxy pairs to- 
gether with details of various physical parameters calculated for 
each system. In Section 3 we present our results and in Section 4 
we discuss these results in the context of earlier work and consider 
the implications for our understanding of interaction-triggered star 
formation. 



2 THE DATA 

2.1 Sloan Digital Sky Survey 

The Sloan Digital Sky Survey (SDSS) is an imaging survey in five 
broad photometric bands (u, g, r, i, and z, defined in Fukugita et alj 
1996) with medium-resolution (R ~ 1800-2100) spectroscopic 
follow up of approximately one million targets. When complete, it 
will cover most of the Northern Galactic Hemisphere. The sample 
used in this work is based on Data Release 1 (DRl) which covers 
1360 deg^ of the sky; however, we use the improved spectroscopic 
data for each object from Data Releas e 2 (DR2). 

The Main Galaxy Sample (MGS. IStrauss et all200 2i) was se- 
lected from the SDSS imaging data as a galaxy sample for spec- 
troscopic observations. Objects were identified as galaxies in the 
SDSS imaging data if their r-band magnitude measured using the 
best-fitting galaxy light curve was at least 0.3 magnitudes brighter 
the magnitude measured using the point spread function model. 
All such galaxies with r-band Petrosian magnitudes brighter than 
17.77 were selected for spectroscopic follow up and were ob- 
served, as far as possible, given the fibre-placement constraints and 
the finite number of fibres on any one plate. The resulting com- 
pleteness is 93%. Of the 7% of galaxies that are missing, 6% are 
galaxies which could not be observed due to the minimum fibre 
separation constraint. 

The MGS was used to construct (Section f2.2t a volume lim- 
ited sub-sample. Each member of this sub-sample was than paired 
with its nearest galaxy (Section f2.4t and its star formation rate es- 
timated using two methods (Section f2.3> . 



2.2 Primary Catalogue 

Our primary catalogue is a complete, volume- and luminosity- 
limited, sub-sample of the MGS. It was de rived from the MG S 
using a procedure analogous to that used bv lG6me~ etal.H2003h . 
The volume of the sample was defined by the redshift constraint 
0.03 < z < 0.1, where the lower limit was used to avoid subse- 
quent large aperture corrections when estimating star-formation 
rates from Ha measurements (Section |2. 3. H . The upper redshift 
limit allows us to construct a complete sample by retaining only 
galaxies with r-band absolute magnitudes Mr < —20.45. 

Two further criteria were needed to remove a small number of 
spurious entries. The first was the elimination of galaxies with red- 
shift confidence less than 0.7. This relatively relaxed constraint re- 
moves most of the galaxies with poorly determined redshifts while 
introducing only a very small bias towards active emission lines 
systems. Secondly, entries in the SDSS with Petrosian z-band mag- 
nitudes fainter than 22.83, i.e., below the detection limit (Table 21 
in Fstoughton et al....2002) . were excluded. Since galaxies in the 
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main sample have r < 17.77, entries with z > 22.83 have extreme 
colours and are almost certainly spurious detections. 

For the purposes of this study, we wish to remove systems in 
which an AGN provides the dominant contribution to the ionising 
radiation (in these cases an accurate estimate of the star-formation 
rate is not possible using th e methods we employ). Following 
IVeilleux & O sterbrockl 1198% . AGN were identified by their posi- 
tions on the [N Il]/Ha vs [O lll]/H/3 and, in a few cases, [O l]/Ha 
vs [O IIl]/Hj8 planes. We required the signal to noise ratio for the 
line fluxes to exceeded 2 before taking them into consideration. 
Approximately 38% of galaxies from the volume limited sample 
were identified as having an AGN. Of these 64% are early- or 
mixed-type galaxies. A further 295 objects which were not clas- 
sified spectroscopically as AGN but exhibited very broad Ha lines 
with CT > 570km s^', were also excluded. 

2.3 Star Formation Rate Estimators 

iKennicutl ^998|) reviews a number of widely used estimates for 
galaxy star-formation rate (SFR), including Ha and [O II] emis- 
sion lines, together with ultraviolet, far-infrared (FIR) and radio 
continuum emissions. Two of these indicators with very different 
properties were used in this work: the Ha line and the FIR contin- 
uum emission. 



2.3.1 The Ha line 

Star formation rates, iff, for all the galaxies in the primary cata- 
logue were estimated from Ha luminosities. Lea, using the rela- 
tion presented in Kennicutt ( 199^: 

wlMr.. yr- 1 ) = . (1) 

'^^^ ' I.27X I034W 

This relation is valid in the case of optically thin recombination, 
and when the star formation obeys the Saloeter ( 1955) initial mass 
function. 

Ha fluxes were calculated from line parameters measured 
by the SDSS spectrosc opic data processing pipeline, outlined in 
IStoughton et alJ i2002h . As discussed above, while our sample is 
drawn from DRI, we used the DR2 spectroscopy data-products 
which are better calibrated and more robust. The pipeline uses a 
number of absorption and emission lines, and broader features, to 
determine redshifts from spectra. Therefore, the Ha line need not 
be prominent to be measured and, indeed, in the majority of ob- 
jects in our volume limited sample it is measured below the noise 
level or in absorption. 

Corrections for dust obscuration, stellar line absorption and 
th e 3" aperture were ap plied following the prescription presented 
in lHopkins et"ai]<2003l) . Both the Ha and H6 lines were corrected 
for stellar absorption. The equivalent width correction s applied 
were 1.3 A and 1.65 A respectively. As discussed in Ho pkins et alj 
i200 3l). these values are lower than the true stellar absorption be- 
cause SDSS spectra resolve the absorption profile. Dust absorption 
correction was calculated from the observed Ha/Hj8 ratio, assum- 
in g an intrinsic ratio o f 2.86 and the extinction curve presented 
in lCardelli et alj ^98^. If the Ha line was measured with a sig- 
nal to noise ratio of less than 2, or if it was measured in absorp- 
tion, no extinction correction was applied. If the measured H/3 flux 
failed these criteria, the extinction was calculated assuming an up- 
per limit for the Hj3 flux of 2.3 x 10^'^ Wm^^ which is approxi- 
mately the detection limit. 

Finally, the Ha fluxes were corrected for the 3"-diameter cir- 
cular aperture of the fibres, considerably smaller than most of the 



galaxies in our sample. The correction was done by scaling the 
Ha flux by the ratio of the total r-band flux to the r-band flux 
measu red within the fibre itself. This procedure has been shown, 
by Hookins et al. (2003), to provide excellent quantitative agree- 
ment with star-formation estimates from global galaxy properties 
such as radio, far-infrared and M-band continuum emission. This 
procedure assumes that the measured Ha equivalent width is rep- 
resentative of the whole galaxy. For the volume-limited sample the 
median physical size probed by the 3" aperture is 4.4 kpc. The fi- 
bre therefore probes a significant fr action of the g alaxy and this 
accounts for the reported accuracy in lHopkins et alJi2bo3) . 

2.3.2 Far-infrared continuum 

Far infrared traces recent star formation, where the natal dust 
clouds still surround the young massive stars and re-processes 
most of their bolometric output. The effects of obscuration at 
these wavelengths are far less prominent than at optical/UV wave- 
lengths, removing a potential source of systematic errors. 

Flux densities at 60 p.m and 100 p.m were obtained for each 
galaxy in our volume limited catalogue by automated query of the 
SCANPI service' , which co-adds individual IRAS scans at the po- 
sition of each of the galaxies. The measurement in each band was 
rejected if the reported signal to noise was less than 2.5 or if the 
template correlation factor was less than 0.9. The FI R luminosity to 
SFR c alibration we applied was the same as used bv lHopkins et alJ 
llOOA namely: 

WiR (Mr^ yr- ' ) = , (2) 

rnKv .J J J 1.39 X 10^6 W 

where 

r 1 + yil^S Lfir > 2. 186 X 10" 
^ I 0-75(1 + Lhr^ 2.186x10^7. 

The effective beam size of IRAS observations is approximately 2 
arcminutes, much larger than the spectroscopic aperture; for this 
reason the FIR data suffer from the opposite problem to that suf- 
fered by the Ha data - namely, that besides accounting for the 
entire FIR output of the primary galaxy, the measured FIR may be 
significantly elevated by the output of other galaxies, most notably 
the identified companion. 

2.3.3 Specific Star Formation Rates 

The star formation rates calculated above are not ideal for study of 
the triggering of star formation as they scale with galaxy mass. A 
more suitable measure is the specific star formation rate, defined 
as the SFR divided by the stellar mass. Approximate stellar masses 
were calculated from z-band (Apff = 9097 A) luminosities, using 
M^(q) = 4.52, as used by Yasuda et al. (2001), and assuming one 
solar luminosity corresponds to one solar mass. The Hubble pa- 
ram eter Ho = 10 km s ~ ' M pc ~ ' is used throughout, for consistency 
with lHopkins et alJ j2003h . 

To calculate the z-band luminosity, we used the reported Pet- 
rosian magnitudes as introduced by Petrosian 1 1976) and discussed 
extensively in the context of SDSS bv lStrauss et aljj2002l) . Study 

' Available at "http://irsa.ipac.caltech.edu/applications/Scanpi/ Provided 
by NASA/ IPAC Infrared Science Archive, which is operated by the Jet 
Propulsion Laboratory, California Institute of Technology, under contract 
with the National Aeronautics and Space Administration. 
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Figure 1. A scatter plot of K- vs. z-band magnitudes for the galaxies in the 
volume-limited sample detected in the 2MASS XSC. 

bv lSchwarzkopf & Dettni^ <200(i) . although of limited size, sug- 
gests that the radial light profile of interacting galaxies is not sig- 
nificantly extended compared to non-interacting galaxies. Even 
if there was a systematic extension, since the fraction of galaxy 
light measured by Petrosian magnitudes depends only on the shape 
of th e light profile of th e galaxy and not on the surface bright- 
ness iStrauss et all20o3) . if the effect of tidal interaction on near- 
infrared galaxy light distribution is to only change its characteristic 
length scale, there should be no bias in estimated masses of inter- 
acting and non-interacting systems. 

The z-band light will inevitably suffer some intrinsic extinc- 
tion. To asses whether it is an un-biased estimator of galaxy mass 
we have correlated our volume-limited sample with the Two Mi- 
cron All Sky Survey (2MASS) Extended Source Catalogue (XCS) 
which is described in Jarrett et al. ( 2000). In FigureQ we compare 
the Petrosian z-band magnitudes with total A'-band magnitudes for 
the subset of galaxies which were detected in the XSC - there is 
a good correlation between them with a dispersion around propor- 
tionality of 0.25 magnitudes. Since A'-band magnitudes are gener- 
ally believed to give an extinction independent estimate of galaxy 
mass we conclude that for this sample the z-band magnitude also 
provides an un-biased estimate. As a further check, we have in- 
vestigated this correlation for close pairs (vp < 50kpc) and find it 
is essentially unchanged from that shown in FigureQwith a mea- 
sured dispersion of 0.27 magnitudes. 

2.4 Defining a Catalogue of Galaxy Pairs 

To investigate the effects of galaxy interactions, we have identified 
the nearest companion galaxy to each member of our primary vol- 
ume limited catalogue. The method used for this was influenced 
by the selection constraints of the spectroscopic follow up sample 
of the SDSS. In particular, the 55" minimum separation between 
targets on any one plate introduces a significant bias against find- 
ing close pairs of galaxies if only galaxies with measured spectra 
are considered. Therefore, we base our search for pairs on all the 
galaxies identified in the imaging data, using measured redshifts 
where available. 

The procedure used was as follows: for each galaxy in the 
primary catalogue, a list of candidate companions out to a pro- 
jected separation of 300 kpc was assembled from the SDSS imag- 
ing galaxy catalogues. Each of these candidates were then, in the 
order of increasing separation, evaluated according to following 
rules until a satisfactory match was found: 

• If the angular separation between the primary galaxy and the 
candidate was less than three times the r-band Petrosian half-light 




2-band magnitude difference 



Figure 2. The distribution (solid line) and cumulative distribution (dashed 
line) of z-band magnitude differences in spectroscopically confirmed pairs. 

radius of the primary (the mean of this value for the volume lim- 
ited sample is 9"), the candidate was rejected. This is necessary to 
remove the small number of cases where the photometric pipeline 
incorrectly de-blends single galaxies into multiple components. 

• If the candidate has a measured redshift, the recessional ve- 
locity difference is evaluated: if it is less than 900 kms^' the can- 
didate is accepted, otherwi se it is rejected. This is a relaxed con- 
straint, which, according to lPatton et alji200oll2002h . includes all 
potentially interacting systems and avoids a bias toward dynami- 
cally bound systems. 

• If no redshift was measured, the z-band magnitude difference 
between the candidate and the primary galaxy is evaluated: if it 
less than 2 magnitudes the candidate was accepted, otherwise it 
was rejected. 

This procedure resulted in 12492 of the 13973 galaxies in the vol- 
ume limited sample being identified with companions. Of these, 
2038 have measured spectra, whilst the remaining systems have 
only imaging data. 

2.4.1 Validity of Spectroscopic-Photometric Pairs 

As described above, candidate pairs were not restricted to the 
MGS, but were taken from all the galaxies identified in the imag- 
ing survey. However, in the absence of redshift information, pho- 
tometric properties must be used to determine if the pairing is just 
a projection effect or if the galaxies are physically close. 

We adopted a relatively simple criterion, comparing the Pet- 
rosian z-band magnitudes of the two galaxies. This technique suf- 
fers both from incompleteness, when there is a large intrinsic lumi- 
nosity difference between true companions, and contamination by 
background galaxies. As both of these effects depend sensitively 
on the cut-off brightness ratio, we now justify the choice of this 
quantity. 

Our choice of Awij = 2 implies that companions which are 
fainter than the primary galaxy by more than about a factor of 6.3 
in the z-band will be rejected, unless they have a measured redshift. 
The z-band was chosen for this comparison as it most closely traces 
the mass of the galaxy. To assess the significance of this incom- 
pleteness we have calculated the distribution of z-band magnitude 
differences in all pairs which have measured spectra, and therefore 
do not suffer from the above effect. The results, presented in Fig- 
ure|2| show that Am, = 2 is in the tail of the distribution and that 
the completeness for this choice of Anij is approximately 80%. 

We have also investigated the probability that a detected pair 
is due to background contamination as a function of the maximum 
magnitude difference. To do this we have used the galaxy count 
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Figure 3. Mean probability of a non-pliysical pair as a function of separa- 
tion from tlie primary galaxy. 

from lYasuda et alj i200lh . The parameterised integrated counts, 
where m denotes the z-band magnitude, are 

A^(z-mag < m) = x iQ^'n-i^) 

deg^ 

The probabihty of a detected pair being due to contamina- 
tion is then calculated as a function of physical separation and av- 
eraged over the whole of the primary catalogue. The results are 
presented in Figure |3| They show that at projected separation of 
over 150 kpc, contamination of our pair sample by background or 
foreground galaxies becomes a significant problem. Therefore, our 
results at separations greater than 150 kpc are significantly diluted 
by non-physical pairs; the reported separations in this regime can, 
however, be regarded as reliable lower limits of separation to the 
true nearest companion. 

2.5 Concentration Index 

Galaxy morphology has long be en known to correlate strongly 
with the star-formation r ate (e.g. ,|Kennic utll|l 998). Additionally, 
numerical simulations bv lMihos & H emguist' 1 1996) show that the 
time dependence of gas inflow during interactions, and the star for- 
mation which results, depend sensitively on the structure of the 
progenitor galaxies, and in particular on the bulge to disk mass ra- 
tio. In order to measure accurately the effect of interactions on star 
formation, it is important to consider these effects. 

The traditional classification of galaxy morphology is the 
Hubble sequence. There is a well established correlation between 
Ha equivalent width and Hubble type, equivalent width increas- 
ing from zero (within observational errors) in E/SO galaxies to 20- 
150 A in late-type spiral and irregular galaxies, with the same trend 
show n between Sa through to Sc spirals (see Kennicutt & Kent 
Il983h . The Hubble classification is based on three characteristics: 
the bulge-to-disk light ratio, tightness of spiral arms and degree of 
resolution of spiral arms. Besides being hard to automate, the Hub- 
ble classification is non-ideal for studies of star formation as the 
last of the three criteria in particular depends significantly on star 
formation. That is, the spiral arms are more likely to be resolved if 
there is significant star formation along them. 

In this work we adopt an alternative morphological classifi- 
cation scheme, the concentrat ion index , which is descri bed in de- 
tail bylMorgai] ilisS) . Okamura et alj ri984) . and Abraham et alj 
il994h . The concentration index, C, is defined as the ratio of Pet- 
rosian 50%- to 90%-light radii as measured in the r-band. Low 
values of this quantity correspond to systems with high central 



concentrations of light which in turn are of an early morpholog- 
ical type. The correlation between classical morphological classi- 
fication and the concentration index f or a sample of ga laxies taken 
from the SDSS is presented in Shima s^u et alj f200ll) . 

We use the concentration index both to exclude early-type 
systems from our study of triggered star-formation and to exam- 
ine the efficiency of this triggering as a function of morphological 
type. For the purpose of removing early-type systems, we exclude 
all galaxies with a concentr ation index less than 0.37 5, larger than 
the 0.33 value proposed bv lShimasaku et all i200ll) . In doing so, 
we reduce conta mination of our late-type sample to less than 5% 
(see Figure 1 1 of IShimasaku et all200lh . at the expense of com- 
pleteness. 

2.6 Pair Morphology 

In the later stages of interaction it is often not possible to resolve 
the nuclei of interacting galaxies, especially in visible light, and 
even when they are resolved the automated algorithm in the SDSS 
photometric pipeline might not de-blend them. More importantly, 
the algorithm used to select nearest companions (Section [2.4> re- 
jects all objects closer than three times the Petrosian half-light ra- 
dius. As a result, our pair sample does not contain the closest pairs 
or merging objects. At the median redshift, and using the median 
half-light radius, the minimum galaxy separation corresponds to 
14 kpc. 

To assess the importance of this effect, we have visually 
inspected three sub-samples, each containing approximately 30 
galaxies. The first consisted of the most actively star-forming ob- 
jects, the second of medium star-forming objects and the third was 
a random control sample. In each case we evaluated the possibility 
that the galaxy is in the process of gravitational interaction via the 
presence of tidal arms and double nuclei. The results are presented 
in Section l33l 



3 RESULTS 

Our primary sample consists of 13973 galaxies of which 12492 
have an identified companion. The distribution of morphological 
types, as determined by the concentration index (Section |2.5> . of 
the the galaxies in the primary sample, together with their com- 
panion, is shown in TableQ The distribution of absolute star for- 
mation rates, Xjf, and specific star formation rate, \jf,„, are shown 
in Figure 13] Here we distinguish between three morphologically- 
defined sub-samples: late-type (C > 0.375), early-type (C ^ 0.33) 
and mixed (0.33 < C ^ 0.375). As expected, given the good corre- 
lation between concentration index classification and Hubble type, 
the vast majority of star-forming systems are classified as late-type. 
The mean star formation rates are 0.7, 1.7, and 3.3 Mq yr"' for the 
early, mixed and late type sub-samples respectively. 

In this and following sections, we consider each pair to con- 
sist of a primary galaxy and companion and examine the star- 
formation rate of the primary galaxy in terms of the properties of 
the interacting pair. 

3.1 Star formation as function of projected separation 

To investigate the dependence of specific star formation rate, \jf„„ 
on true pair separation we consider firstly the dependence of \jf,„ 
on projected separation, looking subsequently at the dependence 
on velocity separation. In Figure |5| we plot the mean specific 
star-formation rate against projected pair separation. The primary 
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Table 1. Distribution across concentration index classes of galaxies in the volume limited sample and their 
companions. 



Primary Galaxy 


Total 


0.375 


Companion Galaxy 
0.33 <C<0.375 


C^O.33 


No companion" 


C > 0.375 


7936 


5120 


1500 


555 


761 


0.33 <C< 0.375 


3344 


2064 


681 


270 


329 


Cs:0.33 


2693 


1455 


576 


271 


391 


Total 


13973 


8639 


2757 


1096 


1481 



" The number of galaxies with no companion identified within the 300 kpc maximum search radius. 



Table 2. The dependence of the specific star formation rate on projected separation 
is characterised by fitting a two power-law function to the specific star formation 
rate vs pair separation data {Fig.|5) after its been divided according to morphologi- 
cal type. 



Concentration 


ai 


/3i 


02 




Index bin 


(10-"yr-i) 




(10-iiyr-') 




0<Cs;0.33 


2.04 


-0.53 


0.11 


1.19 


0.33 < C ^ 0.375 


6.30 


-0.41 


0.70 


0.73 


0.375 <C^ 1.0 


13.6 


-0.04 


2.90 


-1.04 



0.375 < C < 1 
0.33 < C < 0.375 
< C < 0.33 




1-10"' 2- 10"' 3- 

Figure 4. The distribution of (a) absolute and (b) specific star-formation 
rates for galaxies in the volume limited sample, split according the concen- 
tration index into three subsamples. 



galaxy in each pair has been assigned, using the concentration 
index, to one of three morphological classes representing late-, 
mixed- and early-types (see Table0. The data have been binned 
into bins 10 kpc wide for separations less than 100 kpc and bins 
40 kpc wide for separations greater than 100 kpc. For each bin we 
have indicated the mean and the estimated statistical variance of 



the mean; the variance is dominated by sample variance rather 
than measurement errors. In each case we characterise the variation 
of Xjfin with pair separation, rp, by fitting a simple two power-law 
function to the un-binned data: 



Vm{rp) 



■- ai 



25 kpc 



-02 



25 kpc 



(5) 



The best-fitting parameters are detailed in Table|2| 

For each morphological class, there is a significant increase 
in the specific star-formation rate at projected separations less 
than 30 kpc. For late-type systems there is some indication that 
the mean specific star-formation rate decreases systematically with 
projected separation out to separations of about 300 kpc where it 
reaches a value of 1.2 x 10^'" yr^'. Those systems classified as 
early- and mixed-types only show significant enhancement for pro- 
jected separations less than about 25 kpc. The greater dispersion in 
these plots reflects the larger intrinsic range in the estimated star- 
formation rates in these systems, many of which are effectively 
zero. 

The use of the concentration index as a morphological clas- 
sifier does introduce some potential problems. Any enhancement 
in nuclear surface brightness of a galaxy will reduce the value of 
the concentration index for that system. In particular, an intense 
nuclear starburst may result in a late-type system being placed in 
the wrong morphological class. To examine this effect further we 
show in Figure|6|the mean concentration index for all systems as a 
function of projected separation. The concentration index peaks at 
a separation of around 75 kpc and rapidly declines at separations 
smaller than 50 kpc. Inspection by eye of the low concentration 
index systems at small separations, strongly suggest they are disk 
systems with pronounced nuclear starbursts. The enhancement in 
the specific star-formation rate at small physical separations for 
those systems which we have been classified as early- and mixed- 
type is likely, at least in part, to be explained by the misclassifi- 
cation of disc-galaxies with strong nuclear starbursts. We cannot 
eliminate the possibility that there is a real enhancement in star 
formation for intrinsic early-type systems for small physical sep- 
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Figure 5. Specific star formation rates as a function of pair separation for 
late (a), early (b), and mixed (c) galaxy types. 



arations to a companion; further study of this point requires an 
independent means of morphological classification. 

Tidal distortion may also be contributing to the lowering 
of concentration index observed in Figure |6| for close pairs. A 
close encounter between two galaxies can lead to stripping of 
the more peripheral stellar material whilst leaving the central 
stel lar distribution relatively undisturbed, although the study by 
[Schwarzkopf & Dettmar (2000) find no evidence for this. The ef- 
fect of such a distortion would be to increase the value of the Pet- 
rosian 90% light radius whilst leaving the Petrosian 50% light ra- 
dius unchanged, reducing the concentration index. However, ex- 
amination of the relative value of the 50 and 90% light radii with 
pair separation reveal both to decline for close pair separations, in- 
dicating that nuclear starbursts are likely to be the dominant factor 
in the reduction of concentration index at close pair separations. 

The use of Ha luminosity as the star-formation rate estimator 
involves corrections both for extinction and the finite size of the 
aperture as discussed in Section l2.3.1l These procedures may in- 
troduce subtle biases if, for example, the effect of tidal interaction 
is always to trigger a nuclear starburst. To investigate any such ef- 
fects we compare results using the FIR and Ha star-formation rate 



O 0.39 




0.1 0.15 0.2 

Pair Separation (Mpc) 



Figure 6. Mean concentration index, C, as a function of pair separation for 
the full volume limited sample. The dashed line is the least-squares best- 
fitting function of the form C = ai(rp/25kpc)^' + «2(''p/25kpc)'^- ; the 
best fit parameters are: a[ = 1.6, Pi = 0.19, 02 = —1.20, and j32 = 0.24. 



indicators - Figure^ In order to take into account the possibility 
that both the primary galaxy and its companion are contributing to 
the measured fluxes, two mass normalisations of the FIR-derived 
star-formation rates are shown in the figure: normalisation by mass 
of the primary galaxy (as everywhere else); and, normalisation by 
the sum of the masses of primary and companion galaxies if their 
separation is less than 1.5 arcminutes, and just the primary mass 
otherwise (open diamonds on the plot). 

The comparison is limited to those systems for which both 
FIR and Ha are available so this is inevitably biased towards sys- 
tems with higher star formation rates. However, these are precisely 
the systems in which any problems associated with the use of Ha 
as a star-formation rate estimator are likely to be the greatest. We 
find good agreement between the FIR, especially when normalised 
to the sum of the masses of the primary and companion galaxies, 
and Ha data both in terms of the trend with projected pair separa- 
tion and the magnitude of the specific star formation rate. We con- 
clude that the procedure used here to estimate the star-formation 
rate from Ha data gives a good and predominantly unbiased esti- 
mate of the true star-formation rate. 

In order to further investigate the effect aperture correction 
has on our results, in Figure|8|we compare non-aperture corrected 
and aperture corrected star-formation rates as a function of sepa- 
ration; the significance of the aperture correction is made clearer 
by splitting the sample into two redshift bins. It is quite clear that 
the non-aperture corrected data exhibit a similar trend to the cor- 
rected data, the main difference being a somewhat steeper rise of 
the SSFR at close separations. This is consistent with the major- 
ity of the triggered star formation occurring in the nuclear regions 
of the galaxies. Aperture correction brings the two sub-samples 
into good agreement although the higher corrected SSFR of the 
low-redshift subsample in the closest bin suggests that the cor- 
rection may be over-estimated in systems with close companions. 
The implication of this is that galaxies with close companions have 
star formation is enhanced in the nuclear region compared to the 
periphery - that is, they harbour nuclear starbursts. Therefore al- 
though the over-estimate of aperture correction for close pairs con- 
tributes to the uncertainty in our measurement of enhancement of 
star formation in those pairs, the existence of this over-correction 
is in itself strong evidence of nuclear starbursts in galaxies with 
close companions. 

In Figure |9| we show the distribution of specific star- 
formation rates for two sub-samples of late type galaxies: one 
of galaxies with companion within 50kpc projected separation, 
and the other with companions at projected separation greater 
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Figure 7. Specific star-formation rates for galaxies detected by IRAS; open circles: Ha star-formation estimator (as before); solid squares: FIR star-formation 
estimator; open diamonds: FIR star-formation estimator but normalised by mass of both the primary and companion galaxies if they are within 1.5 arcminutes 
of each other 
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Figure 8. Comparison of (a) non-aperture corrected and (b) aperture- 
coiTected star-formation rates versus separation for late-type galaxies in 
two redshift bins: 0.03 < z ^ 0.07 (open circles) and 0.07 <z^0.l (full 
squares). 



than 200 kpc or no companion at all. The Kolmogorov-Smimov 
(K-S) test indicates, at a significance level greater than 99.9%, 
that the star-formation rates in these two sub-samples are drawn 
from different parent distributions. The distributions only differ 
significantly for specific star-formation rates above approximately 
3 X lO^^^yr^' suggesting that the observed increase in the mean 
specific star-formation rate for systems with close companions is 
dominated by systems showing significant increases in their rate 
of star formation. To investigate this further, we show in Fig- 
ure 1101 and Table |3] the variation in Xj/,,, with projected separa- 
tion for late-type systems binned according to the absolute star- 
formation rate: low star-forming (0 < xj/ < SMQyr^'), medium 
star-forming (SMrTjyr^' < iff < lOM^yr^'), and highly star- 
forming (lOM0yr^' < The late-type systems in the lowest 
star formation rate bin < V/ < 3 M0 yr^ ' show no enhancement 
of specific star-formation rate with decreasing pair separation. In 
contrast both of the higher star-formation rate bins (3Mpyr^' < 
xjf < 10M(7)yr^' and 1// > lOM0yr^') show a clear increase in 
specific star-formation rate with decreasing projected separation. 
As is clear from the plots, and the best-fitting parameters given in 
Table |3| this enhancement is most significant at the highest star- 
formation rates. 



In summary, we find strong evidence for an enhancement in 
the specific star-formation rate in late-type systems as a function of 
projected separation to the nearest companion. Within a projected 
separation of 50 kpc there is a very significant increase in i/,„ and 
this can be attributed to the contribution from systems showing the 
largest star-formation rates. 



© 2004 RAS, MNRAS OOO.fTlFni 



Star formation in SDSS close pairs 9 



Table 3. Parameters obtained when fitting a two power-law function to plots of specific 
star formation rate as a function of pair separation for three sub-samples of systems, 
selected according to SFR. 



SFR range (Mgyr"') 


ai(yr ' 


) 


/3i 




) 




0< 


3.26 X 10^ 


11 


-0.10 


2.98 X 10" 


11 


-0.09 


3 < 10 


1.40 X 10" 


10 


-0.10 


1.14 X 10" 


10 


-0.10 


10<\I/^ 1000 


1.43 X 10" 


10 


-1.37 


4.03 X 10" 


10 


0.05 



10"' 



10" 



Tp < SOkpc 
7-p > 200 kpc 
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Figure 9. The relative distribution of specific star-formation rates for galax- 
ies from the late-type volume limited sample with close companions (solid 
line) and no companions or companions which are separated by more than 
200 kpc (dashed line). 
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Figure 10. Specific star formation versus separation for three sub-samples 
selected according to absolute star formation rate: < V' < 3 Mm yr" ' , 
3M0yr"' < ^ < lOM,; yr"' and 1/ > 10M,.,yr"'. 



3.2 Star formation properties as a function of recessional 
velocity separation 

In approximately one-fifth of our volume limited sample the com- 
panion galaxy has a measured SDSS spectrum, and therefore, in 
most cases, a measured redshift. In Figure fTTI we examine the ef- 
fect of the recessional velocity difference on star formation rate 
for these galaxies. The best-fitting parameters from fitting a two 
power-law model are shown in Table |4| The data reveal a decline 
in the specific star-formation rate with increasing recessional ve- 
locity difference out to 900 kms"' . The magnitude of this effect is 
smaller than that observed between specific star-formation rate and 
projected separation. The steepest decline is seen at small velocity 
differences. 

A correlation between projected separation and velocity dif- 
ference may be expected if, for example, all systems were bound. 



2-10-'° - 



iS 1-10- 



200 400 600 800 

Velocity difference (kms-') 

Figure 11. Specific star-formation rate of late type galaxies (C > 0.375) as 
a function of velocity separation. 



Table 4. Parameters obtained when fitting a two power- 
law function to the plot of specific star formation rate 
as a function of the recessional velocity of spectroscopic 
galaxy pairs. 



ai (yr 



j8 



a2 (yr 



4.56 X 10 



-10 



0.05 1.43 x 10" 



-1.38 



We find no correlation between these two quantities for this sam- 
ple; the mean velocity separation is ~ 200km s"' and is indepen- 
dent of projected separation. 



3.3 Optical Morphology 

Optical and near-infrared imaging of IRAS selected galaxies has 
revealed that the fraction of objects which are interacting/merging 
increases systematically with increa sing infrared luminosity. For 
example, I Sanders & Mirabell il99d> find this fraction increases 
from about 10% at log(LFiR/L0) = 10.5-11 to essentially 100% 
atlog(LFlR/i0) > 12. 

As discussed in Section |2!4| our procedure for finding com- 
panions misses very close companions and actively merging sys- 
tems. The results of the previous two sections are subject to this 
bias. In particular, systems which are undergoing strong enhance- 
ments in their star-formation rate due to ongoing mergers are likely 
to appear in our catalogue over the full range of projected separa- 
tions (and/or recessional velocity difference). Such systems will 
therefore only contribute to masking any real effect due to tidal 
interactions. 

We have therefore examined by eye the morphologies of 81 
galaxies in our sample. The 81 galaxies examined were selected in 
three ways: 
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Figure 12. Distribution of specific star nformation rates for galaxies witli 
companions closer than 50 kpc, split according to the morphological type 
of the companion. 

A: The twenty-four most highly star forming galaxies; equivalent 
infrared luminosities log(LFiR/L0) > 11.42. 

B: Twenty-five galaxies with moderately active star forma- 
tion; equivalent infrared luminosities in the range 11.10 < 
log(LFiR/L,3)< 11.13. 

C: A random control sample of thirty-two galaxies. 

Among the 24 galaxies in group A, 38% show evidence of 
being merging systems. This fraction drops to 20% in group B and 
6% in the control sample. The trend observed in the //M5-selected 
samples of increasing optical disturbance with increasing infrared 
luminosity is therefore reflected in our volume-limited sample. The 
large fraction of merger systems at the high star-formation rates 
can only reduce the observed anti-correlation between star forma- 
tion and projected separation. We discuss these effects further in 
Section|3 

3.4 Dependence on the properties of the interaction 

We do not expect all tidal interactions to lead to enhanced star for- 
mation. How effective a given interaction is at triggering gas inflow 
etc. may depend not only on the properties of the perturbed galaxy 
(which we have examined in the previous section), but also upon 
the nature of the perturber. The nature of the perturbing galaxy can 
be probed in two ways: its morphological class and mass. 

In Figure 1121 we show the distribution of specific star- 
formation rates for late-type galaxies with a companion within 
50 kpc. We split the pairs into three sub-samples depending on the 
morphological type of the perturber: late-late, late-mixed and late- 
early pairs. There are no significant differences between these sub- 
samples: the Kolmogorov-Smirnov (K-S) test applied between the 
late-late and late-early distributions gives a formal significance of 
only 20% that they are drawn from different distributions. 

We investigate the effect of companion mass in Figure fTsI 
Again, we split those systems which have companions with pro- 
jected separations less than 50 kpc into three sub-samples accord- 
ing to the z-band magnitude difference between the primary and 
companion galaxy [Am^ = (primary) — m,(companion)]. For all 
specific star-formation rates we find an excess of systems where 
the companion is of lower mass than the primary galaxy. This ef- 
fect is mainly due to form of the luminosity function which rises 
steeply to lower luminosities at the limit of our primary sample. 
Therefore, for real physical pairs, a galaxy of a given luminosity is 
more likely to be interacting with a galaxy of lower luminosity. We 
have used the K-S test to compare the SSFR of the three subsam- 
ples. We tested the null hypothesis that the three distributions are 



Figure 13. Distribution of specific star formation rates for galaxies with 
pairs closer than 50 kpc, split according to z-band magnitude difference. 

all drawn from the same parent distribution. We find no evidence 
to reject the null hypothesis - the formal confidence level that the 
distributions are drawn from different populations is 50%. 



4 DISCUSSION 

The high-quality data from the SDSS has enabled us to construct 
a sample of galaxies with companions which is an order of mag- 
nitude larger than those used in previous studies. Such a data set 
enables us to investigate not only the effect of tidal interactions on 
star formation within galaxies, but also to examine this as a func- 
tion of morphological type and nature of the perturbing galaxy. 
Interpretation of our results is made easier by using a substantially 
complete volume- and luminosity-limited sample. 

We begin our discussion by briefly reviewing sources of pos- 
sible error and bias in our sample. To estimate star-formation rate 
we have used the Ha luminosity from the SDSS. This was cor- 
rected for extinction using measured Hj3 data and then corrected 
for the finite aperture size (see Section l2.3.1> . By constraining the 
lower redshift limit of the s ample to 0.03 t he ape rture correction 
is minimised. As shown in iHonkins et alJ i200.3h this procedure 
leads to measurements of the star-formation rate in agreement with 
other estimators, however, a residual effect remains. For systems in 
which star formation is dominated by a nuclear starburst that falls 
within the SDSS fibre, the aperture correction will tend to overes- 
timate the star-formation rate to some extent. On the other hand, 
if the majority of star formation is off-nuclear, or if the fibre is 
not accurately positioned on the nucleus, then the aperture cor- 
rection will not recover the true star-formation rate. Undoubtedly, 
these effects occur in measurements of individual objects, however 
the data shown in Figure[8|suggest that on the whole the aperture 
correction is a useful way of estimating the global star-formation 
rates. 

We have also made use of IRAS FIR data to give an indepen- 
dent estimate of star-formation rate for those objects in the vol- 
ume limited sample that are detected. These systems are the most 
actively star-forming systems in the sample, and therefore most 
likely to suffer significant problems due to the aperture correction 
and extinction. We have found that the results based on Ha and 
FIR are in good agreement suggesting that the Ha-based estimate 
of the star-formation rate is sufficiently unbiased for the purposes 
of this study. Further interpretation is made difficult by the coarse 
resolution of IRAS; if we normalise the FIR-derived SFR by just 
the mass of the primary galaxy - and ignore the possibility that 
the companion also contributes to the flux - it appears that Ha 
estimated SFR in the close pairs are approximately 30% smaller 
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than those estimated using the FIR data: an indication perhaps that 
we underestimate the extinction in these systems. An argument for 
foll owing this approach is that, at least for samples of bright galax- 
ies. lSurace et alJT2004h and works referenced therein show that in 
pairs of interacting galaxies only one is usually infrared-luminous. 
Alternatively, if we normalise the FIR-derived SFR by the com- 
bined masses of the primary galaxy and companion, we find that 
there is good agreement for the whole range of separations. 

Although we have removed all systems which show any spec- 
troscopic evidence of having an AGN, it is possible that highly 
obscured AGN are present in our sample. Given our procedures, it 
seems likely that if such systems are present their Ha fluxes will be 
dominated by star formation. On the other hand, their FIR fluxes - 
which more closely measure bolometric output ~ may be enhanced 
by the AGN; but, this can only make the correlation between the 
Ha and FIR measured star-formation rates worse. 

As discussed in Section 12.41 closest companions were se- 
lected on the basis of z-band magnitude difference in the projected 
separation range 3x Petrosian half-light radius < rp < 300 kpc, 
and using spectroscopic data in the 20% of cases where it was 
available. This approach has the advantage that our sample is com- 
plete to a minimum pair separation of about 15 kpc and to a well 
defined z-band luminosity difference. This completeness comes 
at the cost of contamination by projection effect pairs; we have 
shown that these only become significant at separations greater 
than 150 kpc. In these cases the distance to the closest companion 
represents a lower-limit to the true separation. Since this contam- 
ination is only significant at larger separations it should not affect 
the interpretation of our results. 

The imposed minimum pair separation and intrinsic difficulty 
in de-blending galaxies with overlapping disks does mean we will 
inevitably miss very close pairs and merging systems. We note that 
of the 9 systems in the highest SFR bin (Section f3.3> that are clas- 
sified as actively merging, all bar one have been erroneously paired 
with distant galaxies. 

We observe an enhancement in the specific star-formation 
rate in late-type systems with a close c ompanion (rp < 50 kpc). 
This resu lt is similar to that found by iBarton et alJ yOOO) and 
Lamba^^L (20Q3) . A quantitative comparison to the results of 
LMnba^ Tal] j20 03|) is possibl e. The functional form used to fit 
the data iii lLambas et alJ i2003h is a single power-law of the form 
(b) = a\r^ where b is the birth parameter which should be pro- 
portional to !//,„. We have fitted our data using a similar functional 
fo rm, restricting the ra nge of projected separations to that used 
bv lLambas et"!?] <2003h [0 < < 100 x (0.7)"' kpc]. After con- 
verting results to a common distance scale we find that we obtain 
agreement with L ambas et al. t2003J only if we apply the addi- 
tional constraint 1//^ > 2 x 10^^' vr~ ', which accounts fo r a bias 
towards more star forming systems in lLambas et alji2003t ). 

The larger sample size, and additional parameters available in 
the SDSS have enabled us to examine these basic results in greater 
detail. 

Although the strongest enhancement in the specific star- 
formation rate occurs for small physical projected separations, we 
also find that for late-type galaxies, there appears to be a system- 
atic decrease of specific star formation with increasing compan- 
ion separation out to 300 kpc. A sim ilar trend may also be present 
in the data of iLambas et all j2003) (see their Figure 1) although 
iT.amhas et alj do not reach this conclusion. This long-range corre- 
lation between star formation and pair separation is seen only in 
late-type systems, i.e., those with concentration indices > 0.375. 

It should be noted that the combination of significant projec- 
tion pair contamination at large separations and the morphological 



separation into only three bins might contrive to make galaxies 
of slightly earlier morphological type biased towards pairs at large 
separations which could produce a correlation similar to that which 
is observed. 

An interpretation of this effect as due to triggering of star- 
formation would, h o wever , be consistent with the predictions of 
Mihos & Hernauist ( 1996), who show that galaxies with weak 
bulge components are more susceptible to formation of bars during 
interaction. The presence of the bar provides the mechanism for 
gas flow towards the nucleus, leading to enhanced star formation 
over a prolonged period of time, even when the galaxies are still 
widely separated. Such weak-bulge systems constitute the highest 
concentration index systems in our sample. In contrast, systems 
with strong bulges (i.e., lower concentration indices) are stabilised 
against bar formation by the presence of the bulge. As a result, the 
onset of star formation is delayed until the final stages of merger. 

It is clear from our data that, if only actively star-forming 
systems are considered (e.g., x^i > lOMoyr^', Figure [To}, the m- 
crease of the average star formation at close separations is much 
higher. This suggests that the significance of interactions will be 
somewhat overestimated in studies which are biased towards se- 
lecting star-forming systems. Our data for the FIR detected sub- 
sample demonstrate this quite clearly. 

We also find a surprisingly tight relationship between the con- 
centration index and pair separation. The steep decline of the con- 
centration index at projected separations below 50 kpc is most eas- 
ily explained if there exists a population of systems with nuclear 
starbursts at these separations. If this is the case, it is quite likely 
that a few systems with strong disk components and intense star 
formation have been misclassified, using the concentration index, 
as intermediate types. It is possible such systems are responsible 
for some, or even most, of the enhancement of SFR in galaxies of 
intermediate type with close pairs. 

The specific star-formation rate is found to decrease with in- 
creasing recessional velocity separation out to the largest velocity 
separations, suggesting that fly-by encounters, in which the two 
systems are not on bound orbits, do contribute to enhanced star 
formation. The specific star-formation rate appears to be a weaker 
function of velocity difference than of projected physical separa- 
tion, however this may be a result of dilution of the velocity results 
due to line of sight projection effects. Comparing our results with 
those of lLam bas et al. ( 200j) we note a similar percentage drop in 
the SSFR over the velocity range 0-200kms^', however, we do 
not observe the marked change in gradient at ~ 350km s^' they 
observe. 

The nature of the perturber does not appear significant in de- 
termining star formation enhancement. We find no difference in 
the effectiveness of triggering star formation between late-type / 
late-type and late-type / early-type interactions. Further, we find 
no evidence for a dependence of the specific star-formation rate 
on the m ass of the companion . This is i n contrast both to th e 
results of 'Donzelli & Pas torizal j 1997b and Lambas et alj j2003b . 
Donzelli & Pastoriza ( 1997j) examine star formation enhancement 
in a sample of 27 physical pairs consisting of a main galaxy and 
a companion that has approximately half the diameter. They find 
that the star formation enhancement is much stronger (~ 50% 
on average) in the minor components of the pairs. In contrast, 
^Lambas et al. (2003), dividing their pair sample into minor and 
major merger candidates according to the relative luminosity of 
the component galaxies, find that the bright components of minor 
merger pairs show a higher probability of having enhanced star 
formation. Comparing the star formation enhancement in major 
mergers, they find no significant difference between the two com- 
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ponen ts. Our results may be reconciled with those of lLambas et alJ 
J2OO3I); as is clear from Figure Elthe mean of the distribution 
for the sub-sample with —2 < Aot. < — 1 is greater than that for 
< Aotj < 2, however this is caused by small number statistics at 
high specific star formation rates. 

An interpre tation consistent with our own results and those 
of lLambas etaP J2003 ) is as follows. The probability of a galaxy 
undergoing enhanced star formation depends strongly on the prop- 
erties of the galaxy itself (availability of gas etc.) and on the nature 
of the orbital interaction with its companion. The mass of the com- 
panion appears to be of less importance - at least up to the factor 
of ~ 6 maximum mass ratio that we probe. There is a strong ten- 
dency for a given system to be involved in an interaction with a 
lower-mass companion (a 'minor' interaction), but we believe this 
results simply from the form of the luminosi ty function. These re- 
sults a re still in conflict with the findings of iDonzelli & Pastorizal 
^993) - in this case it is not clear how their sample selection and 
small number statistics would effect their conclusions. 

In this work, we have not explicitly considered the effect of 
environment on galaxy star-formation rates. Our volume-limited 
sample is selected from the MGS which only contains galaxies 
with spectroscopic data, and therefore is biased slightly against 
galaxies in dense environments, where fibre collisions are more 
likely. Our pair finding procedure however uses all the data, and 
therefore does not suffer from this bias. Thus, although high- 
density regions may be somewhat under-represented in our pri- 
mary sample, the pair separations should be unbiased. 

As noted in the introduction, a number of studies have exam- 
ined the re lationship between star-formation rates and the envi ron- 
ment (e.g.. lG6mez et all2003l : lL'ewis et all20 02: Hashimoto et alj 
Il998h - these studies all find a reduction in the mean star-formation 
rate in dense environments. However, interpreting these results 
is difficult because of the well-known and very strong correla- 
tion between morpholo gical type and environment fe.g.. lDressled 
ll980HGotoetall2003h . a nd the correlation of morphological type 
and star formation (e.g., iKennicutll 1998 1). Other galaxy proper- 
ties are also known to correlate with environment: galaxies are, 
for example, o n average mo re luminous in rich clusters (M^y = 
-20.07 ± 0.07, iDe Propris et al. 2003) an d marginally more lu- 
minous in groups (Mlj = -19.90 ± 0.03, iM^faez^niiSSi 
as co mpared to the field (M^^ = -19.79 ± 0.04. lMadgwick et alj 
l2002h . 

Our results show enhanced specific star-formation rates at 
close pair separations. The simplest interpretation, which we 
favour, is that this is indicative of tidally-triggered star forma- 
tion. Indeed, tidally triggered conversion of gas into stars during 
cluster assembly has been suggested as the cause of the relation 
between density of environment and star formation (Balogh et al. 
I2OO4I and reference therein). Alternatively, it could be a result of 
some other process, also responsible for the observed correlation 
between galaxy properties and environment. This however seems 
unlikely for the following reasons. If close pairs are more likely in 
dense environments, the works cited above suggest that their star- 
formation rate would, on average, be reduced and their luminos- 
ity increased thereby further reducing their specific star-formation 
rate. We, however, observe the opposite effect. Our observations 
could be explained by environmental effects if close pairs are in 
fact more likely in the field where the average star-formation rate 
is higher. A full analysis of these effects is beyond the scope of the 
present paper, but will be the subject of a forthcoming investiga- 
tion. 



5 CONCLUSIONS 

We have used SDSS DRl imaging and DR2 spectroscopy prod- 
ucts to construct a catalogue of nearest companions for a volume 
and luminosity limited sample of galaxies drawn from the Main 
Galaxy Sample. Star formation rates for the volume-limited sam- 
ple have been calculated from extinction and aperture corrected 
Ha luminosities, and IRAS data where available. We used r-band 
concentration indices as a proxy to morphological classification 
and z-band luminosities to estimate galaxy masses. Our results in- 
dicate that: 

(i) The mean specific star formation rate of galaxies shows a 
clear (anti-)correlation with distance to the nearest companion. For 
early- and mixed-type systems this correlation is observed for pro- 
jected separations below 50kpc; for late-type systems, it is ob- 
served out to 300 kpc, the largest separation that we probe. 

(ii) We show this correlation is more pronounced in actively 
star-forming systems. 

(iii) We observe a decline in the mean star-formation rate with 
increasing recessional velocity difference for those galaxy pairs 
where a redshift has been measured for both galaxies. The magni- 
tude of this effect is small when compared with projected separa- 
tion. 

(iv) We observe a tight relationship between the concentration 
index and pair separation which peaks at Vp = 75 kpc and declines 
rapidly at smaller separations. We interpret this trend as due to 
the presence of tidally-triggered nuclear starbursts in systems with 
small projected separation. 
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